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We report temperature- and magnetic field-dependent bulk muon spin rotation measurements in 
a c-axis oriented superconductor CaC6 in the mixed state. Using both a simple second moment 
analysis and the more precise analytical Ginzburg-Landau model, we obtained a field independent 
in-plane magnetic penetration depth A a b(0) = 72(3) nm. The temperature dependencies of the 
normalized muon spin relaxation rate and of the normalized superfluid density result to be identical, 
and both are well represented by the clean limit BCS model with 2A/k J eT' c = 3.6(1), suggesting 
that CaC6 is a fully gapped BCS superconductor in the clean limit regime. 

PACS numbers: 74.70.Wz, 74.25.Ha, 76.75. +i 



Recently the field of graphite intercalated compounds 
(GICs) has regained attention after the discov ery of the 
superconducting GIC CaC 6 with T c ~11.5K.™ Soon 
after the discovery it was suggested 3 that CaC6 is an 
unconventional superconductor, and that the pairing in- 
teraction is of electronic origin. As time elapsed, both 
calculation^ and experiments showed that CaC6 is a 
conventional BCS superconductor. In particular, exper- 
iments revealed the existence of a significant Ca isotope 
effect, 5 the absence of gap nodes P a single gap 7 with 
BCS temperature dependence,^ as well as the BCS tem- 
perature dependence of the London penetration depth P 
Moreov er, ca lculations,^ supported by experimental ev- 
idence l 6 l 7 l 1Q n SD showed the importance of a Ca derived 
interlayer band, crossing the Fermi energy, which have 
sufficiently strong coupling with both in-plane intercalant 
and out-of-plane graphite phonon modes to explain the 
T c of CaC6 within a standard electron-phonon coupling 
mechanism. On the other hand, the coupling between 
graphene-derived electrons and high-frequency graphene- 
derived phonons was also demonstrated to be relevant!^ 

At present the mechanism of superconductivity in 
CaC6 and most of its properties have been investigated 
and clarified. Nevertheless a complete description of the 
superconducting state in a Type II superconductor re- 
quires a whole coherent set of parameters measured with 
high precision, such as the value of the coherence length 
£ and the magnetic penetration depth A. These points 
are somehow still lacking. Indeed, two very different val- 
ues of the in plane penetration depth lambda A a 5(0) have 
been reported in the literatur e 19 * 15 * and the temperature 
dependent superfluid density, p(T), has been measured 
with surface sensitive technique onlyP For the extremely 
air sensitive superconductor CaC6, bulk measurement of 



the superconducting properties would be thus very im- 
portant. We also notice that a new unambiguous bulk 
determination of p(T) could also help in clarifying the 
nat ure of t he CaC6 conduction regime, since both dirty 
limit 9 * 16 * 17 ! and clean limit 15 18 regime have been reported 
so far. This system is also interesting from the point of 
view of /iSR data analysis, since it is a low k = A/£ 
superconductor with reduced, in field, vortex phase. 

The muon spin rotation (/iSR) technique is known to 
be one of the most indicated techniques to measure the 
bulk properties of type II superconductors in the vortex 
state. In this paper we report a transverse field muon 
spin rotation (TF-/iSR) experiment on high quality CaC6 
c-axis oriented samples, aimed to probe the vortex lat- 
tice and determine £ a b(0), A a b(0), and the temperature 
dependence of the superfluid density p(T)J^By using a 
peculiar method developed to collect reliable /iSR data 
on very thin samples, such as the CaC6 ones, we obtained 
A a &(0) = 72(3) nm and (, a b(fy — 38 nm. The penetration 
depth results to be field independent, suggesting a fully 
gapped Fermi surface. The temperature dependence of 
the normalized superfluid density is very similar to that 
one of the normalized muon spin relaxation rate, as it 
is expected for a k ~ 2 superconductor at fields about 
B C 2/2. Both the temperature dependencies are well rep- 
resented by the clean limit BCS model, with 2A/k^T c = 
3.6(1), suggesting that CaC6 is a BCS superconductor in 
the clean limit regime. 

In TF-//SR the sample is cooled-down in a magnetic 
field, which here was parallel to the c-axis. The muons 
are implanted in the sample with their initial spin po- 
larization almost perpendicular to the external magnetic 
field. Each muon precesses around to the local field. 
When the superconducting sample is in the mixed state, 
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each muon probes a slightly different field. As a result, 
a de-phasing process takes place and the average muon 
spin polarization decays. From the average polarization 
of the muon as a function of time, the parameters of the 
vortex lattice, A and £, can be extracted! 20 * 21 1 

CaC6 samples were prepared using the alloy method 
as described by Emery et al.W A stainless steel tube is 
loaded with Lithium and Calcium in the ratio 3:1. Nat- 
ural graphite flakes are then added to the ampoule. The 
reaction takes place in an argon atmosphere for 10 days 
at 350 C. The flakes are typically of a size 2 mm x 2 mm, 
with a thickness of 50 /im. The low-field superconduct- 
ing transition onset is at T c = 11.6 K.^It was reported 
previously that these samples are very reactive and tend 
to degrade very fast. This is observed as an increase in 
the transition width and a decrease in the superconduct- 
ing volume fraction. For that reason the samples were 
kept in an He atmosphere and overall were exposed to 
the room atmosphere for not more than a few minutes. 

The experiments were done in the GPS spectrome- 
ter at the Paul Scherrer Institute PSI Villigen, Switzer- 
land. Three flakes were used, covering together an area of 
about 4 mm x 4 mm. The CaC6 samples are too thin to 
stop the muons. Given the muon stopping distance and 
the density of CaC6, a sample of a few hundred microns 
is needed for this purpose. To maximize the number of 
muons stopping in the samples, we used a method devel- 
oped to measure very thin crystals. 23 The samples were 
sandwiched between two stacks of Kapton foil rectangles. 
This assembly is illustrated in the panel (a) of Fig. [I] 
Each Kapton foil is 125 /im thick. The role of the Kap- 
ton foils in front is to slow down the muons so that most 
of them will stop in the sample. The typical statistics of 
our /iSR measurements was about 17 millions events. 

Muons in Kapton are known to have a very high prob- 
ability of forming muonium (muon-electron bound state) 
However, a muonium signal is not observed in our experi- 
ment. In addition, muons stopped by Kapton do not con- 
tribute substantially to the measured asymmetry since it 
is known that the asymmetry is strongly reduced by a 
thick Kapton layer) 2 ^ The Kapton foils in the back of 
the sample prevent slow muons that did not stop in the 
sample from stopping in the windows of the cryostat. In 
addition, we used a veto detector, that reduces the back- 
ground signal by removing events from the data where 
the incoming muon missed the sample. Using the veto 
counter and Kapton combination we were able to get an 
asymmetry of about 0.1, about half of that we estimated 
to be "good" signal coming from the CaC6 sample. 

The /iSR signal was recorded in the usual time- 
differential way by counting positrons from decaying 
muons as a function of time. The time dependence of 
the positron rate is given by the expression 

N(t)=N —e-^[l + aP(t)]+N hg (1) 
where Nq is the normalization constant, Nb g denotes the 
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FIG. 1: (a) Illustration of the CaC6 sample and the Kapton 
stacks used to degrade the muon beam (see text), (b) Mag- 
netic moment of CaC6 as a function of temperature measured 
at the same field used in TF-/xSR measurements (120 mT). 
The closed (open) symbols represent zero-field cooling (field 
cooling) measurements. Inset to panel (b): magnetic moment 
on enlarged scale to better show the ZFC data. 



time-independent background, r /i =2.19703(4)xl0 -6 s is 
the muon lifetime, a is the maximum decay asymmetry 
for the particular detector telescope and P(t) is the polar- 
ization of the muon ensemble P(t) = j P(B) cos(j^Bt + 
(j))dB. Here 7 M = 2tt x 135.5342 MHz/T is the muon 
gyromagnetic ratio and <j) is the angle between the ini- 
tial muon polarization and the effective symmetry axis 
of a positron detector. P(t) can be linked to the internal 
field distribution P(B) by using the algorithm of Fourier 
transform! 2 -^ 

In Fig. [2] the magnetic field distribution P(B) at an 
applied field of 120 mT above (7 K) and below (3 K) 
T c , obtained from the measured /iSR time spectra by 
performing fast Fourier Transform, is shown. In the nor- 
mal state, P(B) is a sharp symmetric line centered at 
the position of the external magnetic field. Below T c the 
field distribution results to be composed by two signals: 
one is represented by a sharp symmetric peak centered 
at the position of the external field and ascribed to dif- 
ferent sources of background (cryostat window, Kapton 
foils), and the other is a broadened line, which signals 
the formation of a vortex lattice (VL) in the supercon- 
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As a first step, the /iSR time spectra collected below 
T c were fitted by two Gaussian lines 
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FIG. 2: Magnetic field distribution P(B) at an applied field 
of 120 mT above (7 K) and below (3 K) T c obtained from the 
measured /iSR time spectra by means of fast Fourier Trans- 
form. The solid lines are fit by using Eq. [6] Inset: the same 
as in the main panel on larger scales. 



ducting part of the sample. Indeed, VL makes the field 
in the sample highly inhomogeneous. This induces an in- 
crease of the muon spin relaxation rate <r, which is strictly 
related to the second moment of the field distribution. 
Moreover, the line shape of the superconducting P(B) is 
asymmetric, as expected for a field distribution within a 
reasonably well arranged VL. This could be questionable, 
since recently, it has been shown that the irreversibility 
line in CaC6 coincides with the H C 2(T) linep^ which sug- 
gests that the pinning in this sample is very strong and 
that a large number of pinning sites are available. In the 
panel (b) of Fig. [I] the magnetic moment of CaC6, mea- 
sured at the same field used in TF-/iSR measurements, is 
reported as a function of temperature in zero field cool- 
ing (ZFC) and field cooling (FC) conditions. In the FC 
conditions the magnetic moment is negligible compared 
to the ZFC one, indicating a lack of field expulsion. The 
same was observed in the jaSK experiment: we found al- 
most no change in the average muon precession frequency 
when entering the superconducting phase. This is typi- 
cally an effect of strong pinning, which confirms previous 
results P*^ The regularity of the VL could be affected by 
the presence of pinning in the sample. The degree of the 
pinning induced disorder of the flux line lattice (FLL) de- 
pends on many factors, as the strength and the nature of 
the pinning and the value of the external field. However, 
the asymmetric line-shape of the measured field distribu- 
tion in our samples indicates that, although a disorder of 
the FLL is present, it is not large enough to prevent a 
reliable analysis of the data. 
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P(t) = A bg exp(-cr 6 y 2 /2) cos^Bbgt + q 
Aexp(-a 2 t 2 /2) cos(7 /i 5t ■ 

where A^a^, B^g and A,<j,B are the asymmetry, the 
relaxation rate and the mean field of the background 
signal and of the superconducting signal, respectively. 
The first one is assumed to be temperature- and field- 
independent. The second one is an approximation to 
the field distribution in the superconducting state of the 
sample. This distribution is both temperature- and field- 
dependent, and it generates a Gaussian relaxation rate 
<j(T) of the muon spin polarization, which is proportional 
to the second moment of the local field distribution. The 
minimum field for which we could obtain a reliable value 
for <j is 60 mT. Below this field the combination of a very 
fast relaxation and a relatively low frequency prevents a 
good fit. 

In Fig. [3] we show the field dependence of a measured 
at T ~ 2 K. The measurements were done in field-cooled 
conditions: for every field value the sample was warmed 
to a temperature above T c and cooled back to about 2 K. 
As can be seen, there is an almost linear decrease in the 
relaxation rate with increasing field. 

For a fully gapped BCS superconductor, as CaC6 is 
supposed to bej 4 * 6 * 7 * 9 ! the magnetic penetration depth A 
should be field indep endent ! 2 7 | 28 1 In this case, by the anal- 
ysis of the field dependence of a it is possible to obtain an 
estimate of the absolute value of A, by using the approx- 
imation developed by Brandt [Eq.(13) in Refl29]. which 
is considered a very good one for superconductors with 
ft = A/£ > 5 (£ is the coherence length): 



a « 0.1727^(1 
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where 7 M = 2tt x 135.5342 MHz/T is the muon gyro- 
magnetic ratio, $0 is the flux quantum, the in plane 
upper critical field, and A a 5 the in plane magnetic field 
penetration depth. The model reproduces the data rea- 
sonably well (solid line in Fig. [3|. This is thus consistent 
with a field independent A a 5 and, therefore, with CaC6 
being a fully gapped superconductor. The fit (done with- 
out considering the data point at 250 mT which is too 
close or above B C 2 and is compatible with the background 
value) gives A a 5 = 77(3) nm and B C 2 = 230(10)mT. 
The latter, within the Ginzburg-Landau picture, cor- 
responds to £ a & ~ 38 nm. Both the obtained parame- 
ters, £ab an d A Q 5, a re in good agreement with previous 
result J 2 ' 6 ' 7 ' 9 ' 1 ^ ^ and give k = Kb/iab « 2. Although 
the ft = A/£ = 2is rather small the Eq. [3] is still valid in 
the fields ^B/B c2 > 0.5 (see Fig. 6 in Ref 29), that is, 
in the range of our data. 
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FIG. 3: Gaussian relaxation as a function of magnetic field at 
T = 2 K. The solid line is a fit by Eq. [3] 

These results have been obtained by using a gaussian 
approximation for the field distribution P(B). To obtain 
a more accurate value of A a 5 in the following we will 
analyze the muon data by using the analytical Ginzburg- 
Landau (AGL) approximation j 31 l 32 l which is valid for a 
broad range of fields and of Ginzburg-Landau parameter 

K. 

The spatial field distribution B(r) in the vortex lattice 
was calculated as follows: 

B(r) = (B) £exp(-zG • r)B G (A,£, (B)). (4) 

G 

The Bg are the Fourier components! 31 * 32 ^ 

G s (fi + x2 G 2 ) i/2 Kli i tfoo y 

with foo = 1 — 6 4 , and 

& = £(V2 _ 2^)(i + ^4)1/2^ _ 26(1 _ 6)2] i/2_ 

Here, (B) is the average magnetic field inside the super- 
conductor, £ the coherence length, K\(x) is the modified 
Bessel function, r the vector coordinate in a plane per- 
pendicular to the applied field, G = 47r/v'3&(m\/3/2, n+ 
ra/2) the reciprocal lattice vector of the hexagonal FLL 
lattice, a the intervortex distance, and m and n are in- 
teger numbers. From this spatial field distribution the 
probability field distribution P(B) was calculated. 

The /iSR time spectra were fitted with the depolariza- 
tion function P(t) obtained from P(B) as follows: 



cr g describes the disorder of the FLL and was fixed pro- 
portional to 1/A 2 such that cr g (0) = 2.1 /is -1 (note that 
this relaxation is much smaller than the superconducting 
relaxation, which are summed in quadrature). Suc h re- 
lation corresponds to the rigid (well pinned) FLLpHHl 
We neglected additional nuclear relaxation, since the 
concentration of Ca and C isotopes with magnetic nu- 
clear moments is practically zero. The parameter £ was 
kept fixed during the fit at the values determined, at 
each temperature, from the corresponding B C 2(T) via the 
Ginzburg-Landau relation f(T) = vT$o/(27r£ c2 (T))). 
The B C 2(T) values were obtained by linearly interpolat- 
ing the following four experimental points: B C 2(2K) = 
230 mT, given above by the fit of the field dependence 
of a, and B c2 (5.6K) = 120 mT, B c2 (8.5K) = 60 mT, 
B C 2(11.5K) = measured from magnetization experi- 
ments (only the measurement at 120 mT is shown in 
Fig. [I]). These experimental points agree well with the 
B C 2\T) curve reported by Cubitt et al. This model rep- 
resents very accurately the experimental data, as shown 
by the solid line in Fig. [2] 

In Fig. |4j the values of A a 5, obtained by fitting the 
/iSR data with Eq. [6] at the lowest temperature and sev- 
eral fields, are shown. No field dependence of A a 6 is de- 
tected within error-bar, as expected for conventional fully 
gapped BCS superconductors and already suggested by 
the second moment analysis (Fig. [3]). A weighted average 
of the data gives A a & = 72(3) nm, compatible, within the 
errors, with the value obtained by fitting the field depen- 
dence of the relaxation rate, and in very good agreement 
with the value reported in Ref. |9j Recently Cubitt et 
a/P^ reported a magnetic penetration depth A = 50 nm 
using SANS on CaC6 compound. In their analysis they 
used the London model with Gaussian cut off (LGC). As 
it was shown previously by Yaouanc et al., 32 the form 
factors of the LGC model deviate substantially from the 
exact solution of the Ginzburg-Landau model in a very 
broad range of fields. Indeed, an analysis of our data 
by using the LGC model gives A a 6 = 52(4) nm, in good 
agreement with the results of Cubitt et a/P^ The AGL 
model used here gives more reliable results and sys- 
tematic deviations of the AGL model from the exact so- 
lution of the Ginzburg-Landau^ model are within the 
error bar. 

The magnetic penetration depth is defined, in the Lon- 
don limit of low fields, as a measure of the superfluid 
density p oc A -2 . However, for fields close to B C 2 the 
mean magnitude of the Ginzburg-Landau order param- 
eter ip(r) drops due to a substantial overlapping of the 
vortex cores PS The spatial average of the superfluid den- 
sity is thus well described by! 35 * 36 * 



P{t) =A • e-Va^+i* J P(B)e^ Bt dB+ (6) p(T) oc (| ^(r) | 2 )A~ 2 (T) ~ (1 - B / B c2 (T))\- 2 (T), (7) 

Ab g ■ e~ aLtJri<yl ^ BatJr ^ where (• • •) refers to the spatial average over the unit cell 

of the FLL. The temperature dependence of A, needed 
where the second term is the signal of the background to calculate p(T) with Eqj7| was obtained by fitting of 
with asymmetry Abg and relaxation cfl. The parameter Eq. [6] to the data collected at different temperatures in 
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FIG. 4: Magnetic penetration depth as a function of magnetic 
field obtained from the fit to the /iSR time spectra at 2 K. 
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FIG. 5: Normalized superfluid density p(T)/p(fi) and muon 
spin relaxation rate a{T)/a(0) as a function of the reduced 
temperature T/T c measured at 120 mT. Full and dashed lines 
are BCS clean limit and dirty limit curves, respectively. Inset: 
temperature dependence of p and BCS fits. 



a field of 120 mT. The B C 2(T) values were obtained as 
described above. 

In Fig. [5j we show the normalized superfluid density 
p(T)/p(0) as a function of the reduced temperature T/T c 
measured at 120 mT. p(0) is obtained from the BCS fit to 
the temperature dependent p(T) by using the equation!^ 



p(T) = p(0)[l- 



k B T 



poo 

/ f(e,T,A(T))[l-f(e,T,A(T))]de] 
Jo 

(8) 



where /(e,T, A(T)) = [l+exp[^ 2 + A 2 (T) /(J^Tc)]]- 1 
is the Fermi distribution and A(T) = A(0)5(T/T C ), with 
S(T/T C ) given by the conventional BCS temperature de- 
pendence of the gapP^ Indeed, although data at temper- 
atures lower than 1.75 K would be important to clearly 
assess the occurrence of a saturation, however the ab- 



sence of a field dependence of A points to nodeless gap in 
the afr-plane of CaC6, which is consistent with previous 
theoretical end experimental 

studies.™^ Therefore it is 
correct to fit the experimental data using Eq. U In or- 
der to reduce the number of fitting parameters, the fit has 
been performed by keeping T c fixed at the value 5.6(1) 
K obtained from magnetization data at the same field of 
120 mT (see Fig. [TJ. As shown in the inset to Fig. [5| 
the fitting curve (solid line) represents well the experi- 
mental data, yielding p(0) = 85.1(3) pm~ 2 and A(0) = 
0.868(5) meV, which gives 2A(0)/k B T c = 3.6(1). This 
value is compatible with the one expected for a weak 
coupling BCS superc ondu ctor and is in good agreement 
with previous results!^ A fit with the dirty limit BCS 
temperature dependence^ is also acceptable (see inset 
to Fig. [5|, although with larger x 2 and larger errors 
of the fitting parameters (p(0) = 87(2) pm~ 2 , A(0) = 
0.66(2) meV). With these parameters we would obtain 
2A/k#T c = 2.74(1). This quantity is substantially dif- 
ferent from the previously reported one d 6 * 7 * 9 * and it is well 
below the BCS weak-coupling ratio. The normalized su- 
perfluid density, shown in the main panel of Fig.[5j clearly 
follows the clean limit curve, and is not compatible with 
the dirty limit one, suggesting that bulk CaC6 is a BCS 
superconductor in the clean limit regime, in agreement 
with Mialitsin et al\^ and Cubitt et alW^ 

In the main panel of Fig. [5] the temperature depen- 
dence of the normalized relaxation rate cr(T), obtained by 
the two gaussian model (Eq. |2l), is also reported. These 
data closely resemble the p(T) curve. Indeed, as follows 
from the Fig. 6 of Ref. 29 the quantity crn 2 /B c2 (l — b)cx 
(J /[{I — b)/\ 2 ) is practically independent of the reduced 
field b = (B)/B c2 for k ~ 2 and Vb > 0.25. On the other 
hand, cr/[(l — b)/\ 2 ) is proportional to a /p. Therefore, 
the present case with n c± 2 is exceptional, when a sim- 
ple second moment analysis gives a quite precise result 
for the temperature dependence of the superfluid density 
in a broad range of fields. 

To summarize, we performed TF-/iSR measurements 
on the intercalated graphite superconductor CaC6 in the 
vortex-lattice state. Our experimental method included 
the use of Kapton tape to gather good data from very 
thin samples, 50 pm thick. Despite of strong pinning, 
we detected an asymmetric field distribution, typically 
observed in superconducting single crystals in the vortex 
state. The data were analyzed with two different models. 
Using a simple two-Gaussian model, the relaxation rate 
a was obtained. Its field dependence suggested a field 
independent A a fr(0) = 77(3) nm and a coherence length 
£ a b(0) ~ 38 nm. The more precise analytical Ginzburg- 
Landau model yielded A a fr(0) = 72(3) nm, a value re- 
markably close to that one obtained by the second mo- 
ment analysis and in good agreement with previous mi- 
crowave measurements. 9 Also this analysis reveal that 
A a &(0) does not depend on magnetic field, suggesting a 
fully gapped superco nduct ivity in CaC6, in agreement 
with previous reports P 7 -^ The temperature dependence 
of the superfluid density was analyzed by both the clean 
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and the dirty limit BCS models. The experimental data 
are well fitted by the clean limit one, giving a gap-to- 
T c ratio 2A(0) /k^T^ = 3.6(1), in agreement with pre- 
vious results obtained by surface sensitive techniques! 7 -^ 
The temperature dependence of the normalized super- 
fluid density results to be identical to the temperature de- 
pendence of the normalized relaxation rate, as it should 
be for a n c± 2 superconductor at fields B > 0.1 • B C 2. 
Both are well represented by the clean limit BCS curve, 
suggesting that CaC6 is an s-wave BCS superconductor 
in the clean limit regime. 
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